Rashba induced magnetoconductance oscillations in the LaAlOa-SrTiOs 

heterostructure 
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We report measurements of the normal state in-plane magnetoconductance in gated LaAlOa- 
SrTiOa samples. As the orientation of the magnetic field changes within the plane of the interface, 
the signal displays periodic oscillations with respect to the angle between the field and the direc- 
tion of the current. We show that in the underdoped to optimally doped range, a Fermi surface 
reconstruction takes place due to the Rashba spin-orbit term and that the oscillations are due to a 
magnetic field induced opening and closing of a gap at the V point for those Ti out-of-plane orbitals 
having their band minimum close to the Fermi energy of the system. 



In response to a strong worldwide demand for powerful 
portable electronic devices, engineers have designed novel 
architectures for processors while shrinking their sizes 
and slashing down the required energy input. Current 
technology is still mainly silicon-based, but alternative 
strategies are being actively pursued which involve ma- 
terials like graphene or transition metal oxides [iHSI • The 
latter are endowed with an incredible variety of functional 
properties which axe extremely sensitive to structural 
distortions, electronic correlations and crystal chemistry. 
Recently a number of novel phenomena at complex oxide 
interfaces has been discovered |3] and among the driving 
forces that appear to be at play for these new observed 
effects are charge confinement/deconfinement in multiple 
valence heterostructures. Moreover these artificial struc- 
tures can help us control the carrier flow by manipulation 
of the spin. 

One promising candidate is the LaAlOs-SrTiOa 
(LAO/STO) heterostructure. Both compounds are band 
insulators, yet, when a LAO film (4 or more unit cells 
thick) is deposited layer by layer on top of a (001) STO 
substrate, the heterostructure becomes metallic with an 
electron-like character [H [S] . Following the growth pro- 
cedure defined in Ref. [S] one finds that the conducting 
sheet is confined to a few nanometers on the STO side of 
the interface, and that superconductivity sets in below a 
critical temperature T^ in the 0.1 K range [TJIH]- Carrier 
transport was successfully modulated by applying a gate 
voltage V in top- [9 and back-gate geometries reveal- 
ing the existence of a (zero temperature (T)) quantum 
critical point (QCP) for V = Vc separating an insulat- 
ing state from a superconducting state [TU] (see Figure 
1(a)). At low T, in the normal state, charge transport is 
well described by 2D weak localization (WL) in a broad 
range of gate voltages that extends beyond Vc (under- 
doped regime), but it displays anomalous metallic behav- 
ior (AM) for large enough values of V (optimally doped 
regime). Upon increasing V, a steep rise of a Rashba- 



type spin-orbit contribution was evidenced beyond Vc, 
which appears to track the increase in Tc and the onset 
of AM above the dome. 

In this letter we provide further experimental evidence 
of the role of the Rashba spin-orbit coupling on the mag- 
netotransport properties of the LAO/STO interface. We 
present transport measurements for V > Vc, and we 
bring out in full view the remarkable oscillations of the 
magnetoconductance (MC) that take place when a mag- 
netic field B is applied parallel to the interface and then 
rotated within that plane. We correlate the onset of the 
AM with the emergence of a spin-orbit driven Fermi sur- 
face reconstruction. In this regime, as Ti dxz,dyz orbital 
sub-bands are being filled upon increasing V , we con- 
tend that the observed oscillations of the MC are due to 
the periodic opening and closing of gaps at the F point 
of Ti dxz orbital sub-bands having their band minimum 
very close to the Fermi energy. The amplitude of the 
oscillations allows one to extract a spin-orbit energy in 
agreement with the WL analysis. 

LAO/STO interfaces were prepared by pulsed laser de- 
position [11]. A photolithographic process was used to 
define Hall bars for four points DC transport measure- 
ments performed at 1.5 K in a cryostat equipped with a 
rotating sample probe and a 7 T superconducting mag- 
net. Field effect devices were created using the STO sub- 
strate as gate dielectric depositing a gold contact on the 
backside. The current was applied along the [100] direc- 
tion and for perpendicular magnetoconductance the mag- 
netic field was set out of plane along the [001] direction. 
For angular magnetoconductance, the magnetic field was 
applied in plane and rotated from parallel {B ||]100], 
(/)=0) to perpendicular [B ||]010], 0=7r/2) to the current 
(see Figure SI in supplementary material). 

When the applied magnetic field B is perpendicular to 
the interface, the B dependence of the measured MC is 
quantitatively captured using a WL expression for non 
interacting fermions [T] in a broad range of gate voltages 
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FIG. 1: (color online) (a) Evolution of the inelastic (Bi) and 
spin-orbit (Bso) fields as a function of a2o (the sheet con- 
ductance at B=0) for the sample of this work and a reference 
sample (data from Ref. I13[l in the diffusive regime. The top 
voltage scale indicates the corresponding applied voltages for 
the reference sample. The superconducting phase (2D SC, 
blue region) and the transition temperature Tc (triangles) are 
also shown, (b) Evolution of the elastic (Bei) and spin-orbit 
fields with a2D highlighting the crossing at (T2d ~ 0.6 mS. 
Data points extend slightly beyond the diffusive regime, (c) 
variation of the transport effective mass m* with a2D ■ Here 
m*max ~ 2.2me. WL, WAL, AM shaded regions denote the 
weak localization, weak antilocalization and anomalous metal 
regimes respectively. 



extending beyond Vc- The accuracy of the fit allows us 
to estimate the values of the inelastic field, Bi and of the 
spin-orbit field Bgo as a function of the gate voltage. 

Figure 1 (a) displays the evolution of these two quan- 
tities versus <J2d, the sheet conductance (well suited for 
comparing data from different samples). 

Previous studies |13l FH] have shown that the origin of 
the spin-orbit part stems from a Rashba interface term. 
Accordingly, we can determine the effective Drude mass 
m* using the D'yakonov-Perel' relation between the spin- 
orbit (tso) and elastic (r) scattering times ^^ = il^o^T, 
where Ago = hflgo = 2XEkp is the Rashba energy (A is 
the material specific Rashba spin-orbit coefficient, E is 
the interfacial electric field along z, the direction normal 
to the interface, kp is the Fermi momentum) [3] 
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(<I>o = ^ ) and hence r (or B^i the elastic magnetic field) 
in the WL regime (see supplementary material). 



When Eq.ipl holds, i.e. in the range of V such that 
spin-orbit terms contribute to diffusive processes, we con- 
sider that A has little V dependence. By contrast, chang- 
ing the gate voltage modifies the sheet carrier density ns 
and causes a variation of E. The change in E can be 
determined from an electrostatic model of the confined 
electron gas [T^, taking into account the electric field 
dependence of the STO dielectric constant (see Ref. IT7|) . 

Using these prescriptions, we present the variations of 
Bel, Bso and m* with a^jj (or V) in Figures 1 (b) and (c). 
To highlight the trends, we include data points slightly 
beyond the diffusive regime (see below). We note that 
m* is ~ 3.5 lighter in the underdoped regime than in the 
optimally doped regime, which is consistent with the pic- 
ture of an interface charge transport evolving from being 
da^y-dominated (m* ^ 0.64TOe [B], rUe is the bare electron 
mass) to being dxz / dyz-doimnated (to* ~ 2.2me), upon 
increasing V beyond ct^^, - 0.5 mS [131 [Ml [T9H23] ■ The 
concomittant rise in r (Figure S2 supplementary mate- 
rial) fits with this picture, since it signals that additional 
conduction channels start contributing. As reported in 
Refs. (OEDi Tso decreases steeply across Vc, to the ex- 
tent that Tso ~ T for CT°£, ~ 0.6 mS (Figure [l] (b) ) . When 
this occurs spin-orbit processes cease to be diffusive. The 
spin-orbit time is no longer given by the D'yakonov-Perel' 
expression and the Rashba term becomes a bona fide new 
energy scale in the problem, on equal footing with the 
kinetic part. Band structure needs to be recalculated in 
the presence of the (Fermi surface reconstructing) Rashba 
hamiltonian. The evolution towards this AM regime is 
suggestive of a scenario of spin-orbit protected transport 
(against disorder) in the 2D conducting sheet. 

Ab-initio density functional theory (DFT) band struc- 
ture calculations have been carried out for the LAO/STO 
system [51 |23H25]. These studies reveal that the con- 
duction band consists of a manifold of t2g states origi- 
nating from a decomposition of the t2g triplet into sub- 
bands, in response to the interfacial confining electric 
field [SI [T71 127] . The energy minima of these sub-bands 
are located at the F point; the lower the value of the en- 
ergy at the minimum the closer to the interface along z 
the wavefunction is. According to Ref. ([5 ), for values 
of ns ^ 10"'^^ cm~^ the states which lie closest to the in- 
terface, i.e. which do not extend beyond the first 3 unit 
cells from the STO boundary, on the STO side, have d^y 
symmetry and account for a large part of the carrier con- 
centration, dxz, dyz sub-bands spread out deeper into the 
(STO) bulk, and, as they are being filled, they contribute 
to ns on the order of lO""^^ cm~^. These out of plane or- 
bitals lie close the Fermi surface and we argue below that 
they control charge transport in the AM regime. 

To showcase the significant role played by the dxz, dyz 
states, we measure the angular MC of the 2DEG when a 
magnetic field is applied in plane and rotated from par- 
allel to perpendicular to the current. In this geometry 
the magnetotransport is not affected by orbital contribu- 
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FIG. 2: (color online) (a) and (b) Experimental and theoretical plots of the conductance cr2i5(S, <jf>) versus 0, the angle between 
the in-plane magnetic field B and the current for various values of B. (a) corresponds to (j\j^ — 1 mS and (b) to a%j^ — 2 mS. 
(c) and (d) are experimental and theoretical dependence of the conductance on the in-plane magnetic field for (j)=Ti/2. 



tions. Figure [2] shows the variation of the conductance 
with the angle (j) for four values of the magnetic field when 
(T2JJ is 1 mS (a)) and cr2_D is 2 mS (b)). Maxima of the 
conductance are seen for both dopings when the external 
magnetic field is applied perpendicularly to the current. 
However the maximum of the oscillation a2D{B, (p = 7''/2) 
evolves in a non monotonic way for the first doping which 
is not case for the higher conductance state. Note that 
oscillations of the MC in the parallel field geometry have 
been reported by Ben Shalom et al ([55^). For the range 
of sheet conductances and the field intensity that they 
considered (typically higher than ours) they found a pos- 
itive MC for = 0, and suggested that a magnetic order 
forms at the interface. 

For the doping range considered here, we may under- 
stand the behavior of the parallel MC as the field is ro- 
tated within the xy plane if we note that carriers in the 
dxz or dyz sub-bands have one light and one heavy mass 
{nih ~ 20?Tie) [Sllini. Thus, qualitatively, the d^z and 
dyz orbitals are ID-like and, since current fiows along x. 



we focus on the former type which gives a higher con- 
tribution to transport. In the presence of the Rashba 
term, the spin-split bands exhibit an energy gap at the 
F point when B is along x (</> = 0) and a Zeeman-like 
offset when B is along y (^ = f ) (Figure S3 supplemen- 
tary material). The impact of this effect on transport 
depends on the position of Ep relative to Er , the energy 
of the electronic states at the F point when B = 0, since 
the density of states (DOS) at the Fermi energy, g{Ep), 
enters the expression of a2D ■ The conductivity will thus 
show a dip for </) = 0, provided Ep ^ Er ■ If Ep is not 
close to Ey, g{Ep) is almost unchanged as compared to 
its _B = value. As mentioned in the previous section, 
in the underdoped to optimally doped range, Er of most 
of the sub-bands with a d^y character is well below Ep. 
Yet, the Er of dxz symmetry bands are close to Ep. So 
the Rashba induced modulation of g{Ep) is controlled by 
the change in the d^z carrier DOS. Insofar as one may 
reasonably assume that t does not depend on B in the 
2D case, one expects a periodic variation of the MC with 



(f>, displaying crests for B along y and troughs for B along 

X. 

In the framework of this ID picture and with the ad- 
ditional simplifying assumption that the variation of a2D 
with (p is entirely due to the d^z sub-band closest to Ep, 
we find {a2D{B, = tt/2) - ^2^(5, ^ = 0))/a2D{B, cj) = 
0) — l/8{Aso/Ep)'^ when the Zeeman energy is larger 
than the condensation energy of the Rashba state (Ago) 
(see also ^9\). Beyond the qualitative ID model, we have 
modelled the evolution of the band structure for the t2g 
orbitals in the applied B. We use a tight binding model 
featuring kinetic and Rashba terms and we take into ac- 
count the finite value of to/j and the anisotropy of A in the 
xy plane for the d^z, dyz orbitals. The mobility, A and the 
gyromagnetic factor g all depend on V, but we consider 
that they do not change appreciably with nor with the 
magnitude of B in our experiments. For a given V, the 
variation of Sn - i.e. the change in the dxz,dyz carrier 
concentrations - with (j) and B depends on the values of 
the spin-orbit and of the Zeeman energies. The conduc- 
tance <J2d{B,(J)) is then proportional to 5n. Experimen- 
tal data and plots obtained from the model are shown 
in Figures I2] (c) and (d). We note that while a2D{B,(j)) 
increases monotonically with B at fixed (j) for ct^/) — 2 
mS, such is not the case for a2£, = 1 mS, a feature which 
is correctly captured by our model. Figure |3] displays the 
evolution of Aa2D = cr2DiB,(f) = tt/2) — (j2D{B,(f) — 0) 
versus B according to our model and the experimental re- 
sults. We find good agreement using A50 = 7(2.5) meV 
for a 2D = 2(1) mS respectively. These values fall within 
the range of previous experimental estimates |13l I14| . 

Figure |4] shows the evolution of the experimental oscil- 
lation amplitude /S.a2D as a function of U2p, for different 
magnetic fields. As can be seen, A(T2D tends to zero for a 
sheet conductance in the 0.1 — 0.3 mS range. These sheet 
conductance values - which lie in the diffusive regime - 
are close to the QCP, suggesting a potentially important 
role played by the dxz , dyz orbitals in establishing super- 
conductivity. 

In summary, our findings underscore the evolution that 
takes place in the LaAlOa-SrTiOs heterostructure, as one 
tunes the gate voltage in the range where superconduc- 
tivity is observed at low temperature. For low V , con- 
duction is diffusive and is dominated by the dxy orbitals, 
as the impact of disorder is expected to be more severe 
for the ID-like dxzi dyz states. For larger values of V , the 
nature of transport changes and the out of plane dxz , dyz 
orbitals start contributing to the zero field conductivity, 
as evidenced by the evolution of m*, B^i, Bso and by 
the oscillations of the in-plane magnetoconductance. For 
these orbitals, the effect of the strong spin-orbit interac- 
tion has to be taken into account at the band structure 
level for the calculation of their contribution to the trans- 
port. 

We would like to thank M. Lopes and S. C. Miiller 
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FIG. 4: (color online) Evolution of the amplitude of the ex- 
perimental oscillations Aa2D {Aa2D ~ 0-20 {B,<l) = f ) — 
<72d{B, (j) = 0)) as a function of (T213 ■ "Plus" symbols pertain 
to measurements performed on an additional sample at 7 T. 
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SAMPLE PREPARATION AND EXPERIMENTAL GEOMETRY 

A LaAlOs (LAO) layer thicker than 4 unit cells was grown on a (001) Ti02-terminated SrTiOs (STO) substrate 
heated at 800°C in an oxygen pressure of lO"** mbar. The KrF excimer laser fluency was set to 0.6 J/cni^ with a 
repetition rate of 1 Hz. After deposition, the oxygen pressure was raised to 0.2 bar and the sample was kept at a 
temperature of about 530°C during 1 hour before cooling to room temperature. The layer thickness was estimated 
in situ from intensity oscillations of the reflection high energy electron diffraction (RHEED) spots. 




FIG. 1: Shematic view of a field efi'ect device. is the angle between the current (J) and the in-plane magnetic field Bin-piam- 



WEAK LOCALIZATION ANALYSIS 

In the diffusive regime, the perpendicular magnetoconducance (MC) is quantitatively captured using weak localiza- 
tion corrections. We have estimated the spin-orbit fleld Bgo and the inelastic fleld _B,; as a function of the gate voltage 
using the Maekawa-Fukuyama expression for the magnetoconductance [T]. In order to obtain both the effective mass 
m* and the bare (Drude) elastic transport time r, we need to determine the bare conductance (ct^) for various V. 
To that end, we use the experimentally measured B = sheet conductance ((721)) ^^'^ substract the zero field WL 
correction to the conductance, given by 
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In order to determine to*, we consider the D'yakonov-Perel' relation |2j between Tg^ and r, namely ^=^ — i^l^r 
where Ago = fi^so — ^XEkp is the Rashba energy [E is the interfacial electric field, kp is the Fermi momentum). 
In the diffusive regime, the spin-orbit fleld is deflned as Bg^ — -r^ — , where <&o = y^ and D — hvp^r is the Drude 



f and L> : 
(bare) diffusion constant. Using the D 'yakonov-Perel' relation, it follows that Dtso — ^^_F^jf|-- Taking into account 



the fact that nfvp — fikp and that Qso — '^^^^'' we obtain 
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Elastic scattering time 



The elastic field in Eq.J2| is given by i?e( = 



4^£)-. Using D — ^vp'^T, m*vp = hkp, and kp = lixn allows us 
to rewrite B^i in terms of the Drude mobility /^jj = -^ and of the mobile sheet carrier density n. Because Acr^^/^ 
is a correction to the conductivity and because B^i enters Eq.([2| through a Log term, we may replace /i£) by the 
renormalized /x - obtained from measuring the sheet conductivity - in the expression of Be; when we compute Acr^^ . 
We then subtract Acr^^ from the measured sheet conductance o'^y) to get the Drude conductivity ct^ = ne^D. 
Besides, magnetotransport allows us to extract the sheet carrier density and the effective mass (Eq.([3|); from these 
three quantities we determine the elastic scattering time r and its variation with (72£) is shown in Figure v2\ We note 
the steep rise in r which mirrors the sharp increase in m* beyond a2D ~ 0.5 mS. Such upturn may be understood 
if we recall that while d^y orbitals appear to dominate transport for lower values of <72d- d,xz,yz states contribute 
significantly for larger (T2j3.The steep rise suggests that the mobility of the dxz.yz orbitals is sizably larger than that 
of the dxy orbitals. 
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FIG. 2: Dependence of the elastic scattering time (left scale) and of the effective mass (right scale) on the sheet conductance 



BAND STRUCTURE MODEL FOR IN-PLANE FIELD MAGNETOCONDUCTANCE 



Beyond a characteristic voltage V > Vc, t^o can become less than r - the transport elastic time - implying that the 
spin-orbit term is no longer diffusive but that it changes the eigenenergies of the states. Figure |3] displays the band 
structure for the d^z orbitals split by the Rashba coupling. 




FIG. 3: Shematic view of the band structure according to our model, a) We show a d^z sub-band spht by a spin-orbit couphng 
of 5 10~^^ eV m. The effective mass used in the calculation is m*^2=0.64 me. -Br defines the crossing point of the spin-split 
bands, b) Applying a magnetic field parallel to the current direction ([fOO] direction, = 0) opens a gap at the gamma point: 
when the Fermi level is at Er, the change in carrier density leads to a decrease of conductance, c) When the magnetic field is 
in plane and perpendicular to the current ([010] direction, (j) — 7r/2), a Zeeman effect is observed on the band structure. 



MODELLING PARALLEL-FIELD MAGNETOCONDUCTANCE OSCILLATIONS 

In order to describe the evolution of the 2DEG conduction bands as the orientation of the parallel field B changes 
in the plane of the interface, we use a tight-binding representation of the t2g triplet states. In nionientuni space, the 
Haniiltonian has diagonal entries H^ consisting of 2 x 2 spin blocks; H^ contains a kinetic part 



(2i^(l - cos (k^d)) + 2ty{l - cos (fcyd)))H^2 



(4) 



a Rashba piece 
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and a Zeeman term 



±> COS0 Ox i> sin0 ay 
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H^2 is the unit matrix in spin space and cr^,, cTj,, cTz are the Pauli matrices. The angle between the in-plane magnetic 
field B and the x (current) direction is (p. d (~ 3.91A) is the lattice spacing. 

dxy orbitals are invariant under a, x ^ y symmetry, so that tx = ty and axR = ayR = aR. By contrast, dxz and 
dyz states are anisotropic and, from the small k limit, ty/tx (tx/ty) is equal to the ratio of the light mass mi over 
the heavy mass mh for the dxz (dyz) orbitals respectively. Similarly, the coefficients of the Rashba contributions are 
anisotropic and, following the discussion presented in Ref. [3,, we account for this by setting UyRJaxR = mi/mh 
{axR/oiyR — mi/rriii) for the dxz (dyz) orbitals respectively. 
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To the above energy we must add the energy term describing the confinement in the direction perpendicular to the 
interface. It produces an a priori different value of the energy Er at the F point for each of the d states [B]. With 
these energies, we may determine the occupations of the orbitals for a given value of Ep- When current flows in 
the X direction, we focus on the dxz orbital whose -Br is close to Ep. If the elastic transport time does not change 
appreciably upon varying B and (/), we claim that the angular dependence of the magnetoconductance stems from the 
modulation of the DOS as B rotates within the plane of the interface. So we write a2D{B,(j)) = a'^^\n) + C{n)Sn. 



a^^^ (n) and C{n) depend on the mobile carrier concentration n (on the gate voltage) but not on B or (j). The carrier 
density in the dxz band whose Er is close to Ep, 5n (<< n), changes with the gate voltage, B and (f>. 

a^^\n) is the contribution to the conductance of those d states whose Er is not close to Ep- As explained in the 
text, we expect negligible B and (j) dependence in that case. 

For a given ct^j^ we choose one value of the field (7T for (T2d ^ 2niS, 5.5T for ct^d = IniS) and we determine 
a^^'{n), C{n), a^j and g so as to best fit the experimental variation of (72d{B,4>) with (f>. We subsequently use these 
same values when we change B and compute the new a2D{B,(j)). A comparison between the theoretical and the 
experimental curves is shown in Figures 2 and 3. 
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